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Abstract 

By using a small rotation approach, we show that it is possible to 
obtain well behaved perturbed solutions for the amplitude of the electro- 
magnetic field propagating in a photonic waveguide array. 

1 Introduction 

The analogy between linear lattices and atom-field interactions PQ has been a 
fundamental step for the emulation, via classical interactions, of quantum me- 
chanical systems. This is not only important because of pure scientific reasons 
but also because of the possible applications in quantum information processing. 
In this latter case, the properties of classical systems have been used to realize 
quantum computational operations by quantum-like systems and, in particular, 
it has been show how a controlled-NOT gate may be generated in nonhomo- 
geneous optical fibers [2J. At a fundamental level, e.g., it has been possible to 
the emulate the most basic atom-field interaction, the Jaynes-Cummings model, 
theoretically |3] and experimentally [4] with arrays of photonic waveguides; and, 
just to give another example, it has been proposed to model non-linear coherent 
states [5] in waveguide arrays [6] . Linear coherent states have also been modeled 
via linear arrays of photonic waveguides [?]• 

2 Diatomic waveguide arrays 

In a diatomic photonic waveguide the field amplitude in each way is obtained 
form the infinite system of coupled ordinary first order differential equations [3] 

du 

i—r^ - w(-l) n u„ - a(u n+1 + u n -i) — 0, n= -co, ...,oo, (1) 



where a and u) are arbitrary constants. The field at z — may be given, in 
general by u n (z = 0) = tp n , for which we only ask to be normalized. 
One can associate with this system a Schrodinger-like equation 

i—^—=H\i/j(z)), (2) 

subject to the initial condition \tp(z = 0)) = Y2m=-oo tpm\ m )- The "Hamilto- 
nian" of this Schrodinger-like equation is given by 

H = uj{-lf + a{V + V*), (3) 

and the linear operators V and V* are defined as 

oo oo 

V= \n)(n + l\, Vt= J2 \n + l)(n\, (4) 



n— — oo n— — oo 



acting over the vector space generated by the complete and orthonormal set 
{\n);n — — oo, +00} . 

If we propose as solution of the Schrodinger-like equation ([2]) the expansion 

00 

|#0>= E u n(z)\n), (5) 

k— — 00 

we obtain for the coefficients u n (z) — (n\ip(z)) the infinite system of first order 
ordinary differential equations (p}. Thus, by solving the Schrodinger-like equa- 
tion [2] means to solve the system [T] 

The action of the operators V and on the states \n) is 

V\n) = \n-1) and F t |n) = |n+1). (6) 



3 Exact solution 

The formal solution to the Schrodinger-like equation is 

\^) = e - iz6 \m), (7) 

where the Hamiltonian H is given by ([3]) and where we have already used the 
initial condition IV-'(O)) = \ m )- 

Because of the commutation relations of the operators involved in the Hamil- 
tonian is not simple, one can not give a closed form for the evolution operator. 
In order to give a solution, we will make a change of basis, we will go from 
the discrete basis {\n); n — —00, ...,00} to a continuous basis denned by (the 
Fourier series) 

00 

\4>)= J2 e m V>- (8) 



2 



The inverse transformation is clearly 



\n) = ±- I #e-^|^). (9) 



2tt 

Therefore, the solution to the Schrodinger-like equation can be written as 

\1>) = r d^e-^e^l^). (10) 

27r J -TV 

Thus, we have to analyze only the last part of the above equation, the action of 
the operator e~ ltH on \<f>). We write the definition of the exponential operator, 
and after we split the even and odd powers, 



k. — A: — i \ 

fc=0 v y fc=0 v ; 

We have to study now the action of the Hamiltonian H and the square of the 
Hamiltonian H 2 on the wave functions \<p). Using the definition the fact 
that the operators V and commute, and that they anticommute with the 
operator (— 1)™, it is easy to show that 

H\4>) = w|0 + 7r) + 2acos0|0), (12) 

H 2 \(j)) = ft 2 |0), (13) 

and 

ff 2 |0 + 7r) = f> 2 |0 + ^), (14) 

where 

il(4>) = yjb) 2 + 4a 2 cos 2 (j>. (15) 

With all these results, we get back to equation (fTT|) . to obtain 

e-" 6 \<f>) = cos(ttz)|0) - i 2acos0^j^|0) - ioj^^-^ + tt). (16) 



Using now expression (|10p . we get 



.. . , sin(51z) sinffte) . , . 
COS(S22:)|0) — i2aCOS0 — — — \<p) — lb) — 10 + 71") 



To obtain the solution to the infinite system of differential equations, we recall 
that by definition u n = (nltp), so 

u n (z) = ^f d<Pe~ lm * (17) 

f r- / i \ i / i ; \ „ , sin[f2((/))zl , . ,. sin[f2((/))zl , . , . 1 

x |cos[Q(0)z](nl0)-i2acos0 qAV < w ^ + 7r >) 
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It is very easy to prove that (n\<p) — e m *, and substituting this in the above 
expression, we get the final answer to our problem, 

u n {z) = ^f dc^e-^-^ (18) 
f r^/ ,•, i „ s'm\£l(6)z} , .„ s'm\ft(d>)z\ 

X |ooB[n(^]-^OOB^-^-(-l)- W -^ 

As the waveguide array is symmetric and infinite, we do not loose generality at 
all if we consider that the " central" waveguide is shined (m = 0) . If the initial 
state is \m) — |0), the result is reduced to 

( \ 1 f .«J \cmj.\ i n sm{n((j))z} sin[O(0)z 

u n (z) = -J ^ |cos[Q(^] - *2acos^-±^-i - (-iJ^-I-L^ 

(19) 

Using the parity properties of the functions involved, we can finally write the 
solution of our original infinite system (|T|), when m = 0, as 

u n {z) = - [ cosincj)) f oos[fi(0)«] - i [2acos^+ (-l) n w] sm ^^ z ] 



n J { 1 fi(<£) 

(20) 

Note that this solution satisfies the initial conditions. Indeed, if z = we get 
u n (z) = ^ J cos(n(p)d(j) that it is equal to 5 n fi- Also note that if there is no 
interaction between the guides; i.e., a = 0, we obtain 

u n (z) = e-^-^S^o (21) 

that it is the solution of the trivial system of equations obtained. 
In Figure 1, we present the behaviour of the exact solution of the guides to 
10, with lo = 1 and a — 0.3 for z from to 100. For symmetry reasons the 
behaviour of the guides —1 to —10 is exactly the same. 



4 Small rotation 

We consider the case in which the parameters obey a <C w, and transform the 
Hamiltonian 

H = u(-lf + a(V + V*), (22) 
via the unitary transformation 

R = e^- 1 ^^ + ^\ flt = e -fe(-i) ft (V+V t ) j (23) 

such that H R = RHB). We use the relation e sA B e - sA = B + s[A,B] + 
^r[A, [A, B]] + . . . and note that we can cut the series to second order because 
a <C u> [5] , obtaining 

H R « u(-lf + ^(-lf(V + V*) 2 . (24) 
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100 * 10 



Figure 1: Exact solution of the guides to 10 (for symmetry reasons the behavior 
of the guides —1 to —10 is exactly the same) with uj = 1 and a = 0.3 for z from 
to 100 

Because the parity operator (—1)™ and the operator (V + V*) 2 commute, we 
can obtain the propagator operator, Ur(z) — exp(-iHRz), as the product 



(25) 



and the propagator associated to ([22]) is then 

u(z) * £t e — z(-if e ^4f(-if(v+v<r^ 



(26) 



Expanding the square in the exponential, using once more the commutativity 
between (—1)™ and the squared V operators and after some trivial algebra, we 
can write the ket \fp(z)) as 



\<ip(z)) = i?V 



(-1)' 



e 2 " 



(27) 



with |^(0)) the initial condition related with the waveguide that is shined. As 

- r . » -| -1 

i(— l) n (V^) 2 = —i(—l) n V 2 , we may develop the second exponential above 

in terms of Bessel functions by using their generating function and obtain 



k=—oo 



)(-!)" 



[(-if] k k\m) 



(28) 

(—1)"^ , again the properties of the V operators 
and the generating function of the Bessel functions, we can write the following 



Using that (-l) fl V r t 
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Figure 2: Comparison between the exact numerical solution (continuous black 
line) and the small rotation approximation solution (dashed red line) for u = 1 
and a — 0.1, and for the 3 first guides 




Figure 3: Comparison between the exact numerical solution (continuous black 
line) and the small rotation approximation solution (dashed red line) for w — 1 
and a = 0.3, and for the 3 first guides 



explicit expressions for the operators R and W , 

oo 



(29) 



j = -oo 



G 



and 



fJi— — oo 



(30) 



that after being substituted in Equation (f28j) give us 



oo oo oo 



W*)> = E E E H) fc (-i)^ 

k— — oo j — — oo fl — ^OD 



2 

a Z 



J, 



(31) 



v 2k \ i (-iy i v\ < 



[{-iff (-iff |^(0)). 



It is possible to show that 

(-If V " |m) = (-1) 



- jm— ^ J ' 



\m - j) 



(32) 



for j positive and negative; so considering the initial condition \ift(0)) = \m), 
and after some algebra 



u n (z) = {n\ip(z)) m = (-1)' 



■»(m-l)- n (n-l) 



OO OO , 

E E (-l) fc(m ^e _ ' (_ir "T 



fc— — oo j= — oo 



... , / a 2 z 
i k J } 



Jj ( — J Jn-m+2k+j 



is obtained. 

If we consider that the "central" waveguide is shined (m = 0), we get 



u n (z) = (-1)" ™ 2 1 



OO OO 

£ 2 (-l)*e- <-'»'( 



k—-~ oo j — — oo 



that is the final solution in this approximation. 
Using that 



■Jn+2k+j 

(33) 



E^(z)] 2 = kl-[^o(z)] 2 } 



(34) 



and that 

oo 

E JiWA+iW = 0, (35) 

j=— oo 

it is possible to show that the initial conditions are satisfied. It is also easy to 
verify that in the special case when a = 0, we get the correct very well known 
trivial solution. 

In the Figure 2, we compare the exact solution and the small rotation approxi- 
mation for a = 0.1 and in Figure 3 for a — 0.3 giving still good accuracy even 
for larger a's. 
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5 Conclusions 



It has been shown how we can obtain, to a good degree of accuracy approx- 
imated solutions to the problem of light propagating in waveguide arrays. A 
Rayleigh-Schrodinger perturbation approach had been proposed previously by 
Rother, [TU] that however gives rise to convergence problems. Here because of 
the method applied, i.e. the small rotation approach, these convergency prob- 
lems are eliminated. 
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